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ABSTRACT 

We have observed the galaxy environments around a sample of 21 radio-loud, steep- 
spectrum quasars at 0.5 < z < 0.82, spanning several orders of magnitude in radio 
luminosity. The observations also include background control fields used to obtain the 
excess number of galaxies in each quasar field. The galaxy excess was quantified using 
the spatial galaxy-quasar correlation amplitude, B m , and an Abell-type measurement, 
A0.5 (Hill & Lilly 1991). A few quasars are found in relatively rich clusters, but on 
average, they seem to prefer galaxy groups or clusters of approximately Abell class 0. 
We have combined our sample with literature samples extending down to z 0.2 and 
covering the same range in radio luminosity. By using Spearman statistic to disentan- 
gle redshift and luminosity dependences, we detect a weak, but significant, positive 
correlation between the richness of the quasar environment and the quasar's radio 
luminosity. However, we do not find any epoch dependence in B m , as has previously 
been reported for radio quasars and galaxies. We discuss the radio luminosity-cluster 
richness link and possible explanations for the weak correlation that is seen. 
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1 INTRODUCTION 

Since the discovery of quasars, much effort has been ex- 
pended in trying to understand the physics of active galactic 
nuclei (AGN). Most work has concentrated on their nuclei, 
but during the last 10-15 years the galaxy environments 
around different classes of AGN have received much atten- 
tion. 

The first studies that were made in this field (e.g. 
Stockton 1978; Longair & Seldner 1979; Yee & Green 1984; 
Yee & Green 1987 (hereafter YG87); Prestage & Peacock 
1988, 1989; Yates, Miller & Peacock 1989) showed that ra- 
dio quasars and galaxies are located in regions of enhanced 
galaxy density. These findings were confirmed by more re- 
cent studies (e.g. Ellingson, Yee & Green 1991 (hereafter 
EYG91); Hill & Lilly 1991; Allington-Smith et al. 1993; 
Fisher, Bahcall & Kirhakos 1996; Zirbel 1997), which added 
more weight to the consensus that the very local environ- 
ment may be important for the formation and evolution 
of AGN. It is believed that AGN formation is closely re- 
lated to the formation of supermassive black holes, possibly 
by galaxy-galaxy interactions and mergers. Interactions and 
mergers may also provide the fuel supply to a pre-existing 
central black hole. The study of AGN environments may 



therefore provide important clues about their formation and 
evolution. 

Radio-loud AGN are found to be associated almost ex- 
clusively with giant elliptical host galaxies (e.g. Taylor et 
al. 1996; Lacy et al. 1999 and references therein). Both the 
need for interaction, and the apparent need for a massive 
host, make it not surprising that groups or clusters of galax- 
ies might be associated with radio-loud AGN. Furthermore, 
because the efficiency of conversion of the bulk kinetic en- 
ergy of the radio jets to radio luminosity increases with en- 
vironmental gas density, selection effects will also tend to 
favour finding steep-spectrum radio sources in gas-rich envi- 
ronments, i.e. in groups or clusters. High-jz radio-loud AGN 
can therefore be used as tracers of groups and clusters at 
early epochs, and the evolution of galaxy groups and clus- 
ters can thus be studied independent of the selection biases 
associated with optical or X-ray cluster surveys. 

Studying AGN environments may also help establish 
the relationship between different classes of AGN. In partic- 
ular it has been proposed that radio galaxies and radio-loud 
quasars are the same type of object, but viewed at different 
orientations to the line of sight (Scheuer 1987; Barthel 1989; 
Antonucci 1993). The most robust test of such orientation- 
based 'unified models', is to select the AGN in question with 
one orientation-independent quantity (e.g. radio lobe lumi- 
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nosity) and compare another (e.g. environment). And in- 
deed, observations of fields around radio- loud quasars and 
galaxies show that they are sited in environments that oc- 
cupy approximately the same range in galaxy density (e.g. 
Longair & Seldner 1979; Smith & Heckman 1990; Hill & 
Lilly 1991; Yee & Ellingson 1993). 

For obvious reasons, most studies have been concerned 
with sources at 2<$.6, but recently various lines of evidence 
point to radio-loud AGN being in cluster environments also 
at *£1 (e.g. Hall & Green 1998; Dickinson 1996). Compli- 
cations due to the possibility of lensing bias (Hammer & 
Le Fevre 1990; Fort et al. 1996; Benftez, Martinez-Gonzalez 
& Martin- Mirones 1997; Schneider et al. 1998), the possibil- 
ity of non-thermal contributions to X-ray emission, and the 
low contrast of clusters against the background, however, 
make detections of clusters at these redshifts hard to prove. 

An apparent decrease in the galaxy density around pow- 
erful (FRII) quasars from z ~ 0.5 to the present epoch has 
been observed (YG87, EYG91), although the amount of evo- 
lution is still not firmly established. Several investigators 
have also found that powerful FRII radio galaxies inhabit 
richer galaxy clusters at z ~ 0.5 than at z<,0.3. On the other 
hand, the less powerful FRI galaxies seem to be commonly 
found in clusters at both higher and lower redshifts (Lon- 
gair & Seldner 1979; Prestage & Peacoc k 198 8, 1989; Hill 
& Lilly 1991; Zirbel 1997). Yates et al. (jl98§) studied the 
environments of powerful radio galaxies in a wide redshift 
range and found th at sources at z > 0.3 reside in clusters 



of typically Abell (Abell 1958) richness class 0, while the 
enviro nments at z < 0.3 are three times poorer. Hill & Lilly 
(1991) compared a sample of radio galaxies at z ~ 0.5 span- 
ning a wide range in radio luminosity with a lower redshift 
sample by Prestage & Peacock (1988, 1989) and confi rmed 
the Yates et al. result. Allington-Smith et al. (1993) also 
arrived at a similar result, but with a less strong correlation 
between redshift and richness. 

Most of the samples that have been studied are flux- 
limited, i.e. they suffer from a strong correlation between 
redshift and radio luminosity. This raises the possibility that 
the observed evolution in richness is dependent on radio lu- 
minosity rather than redshift. Yates et al. (1989) believed 
that this was the case in their study, but they were unable to 
confirm it since their sample had this strong redshift-radio 
lumin osity correlation. However, for radio galaxies, Hill & 
Lilly (1991) were able to show that this is not the case by 
observing z ~ 0.5 sources over a wide range in radio lu- 
minosity. They found no correlation between richness and 
radio luminosity within their sample. 

In this paper, we investigate the situation at 0.5 < z < 
0.82 for radio-loud quasars. This redshift range was cho- 
sen to extend previous work, which went up to z ~ 0.6, to 
as high a redshift as possible consistent with keeping the 
redshifted 4000 A break shortward of the /-band. We se- 
lected quasars randomly from flux-limited samples spanning 
as wide a range in radio luminosity as possible. By extending 
the redshift range of previous studies, whilst maintaining the 
luminosity range, we are able to better disentangle trends in 
environmental richness due to cosmic evolution from those 
due to radio luminosity. The work presented here is also part 
of a larger study which aims at investigating and comparing 
the galaxy environments of different classes of AGN within a 



relatively small redshift range, but over as wide a luminosity 
range as possible. 

We quantify the amount of galaxy clustering in each 
quasar field using the spatial quasar-galaxy cross-correlation 
amplitude, B gq , and an Abell-type measurement, iVo.5- 
These two quantities are briefly described in Section ^. In 
Section |3| and ||, we present the quasar sample and describe 
the observations. The data reduction, photometry and faint- 
object detection in the CCD images are discussed Section ^| 
In Section ^, we quantify the galaxy environments in the 
quasar fields using _B gq and iVo.s, and compare our results 
with previous work on radio-loud AGN and other two-point 
correlation functions. In Section ^, we investigate redshift 
and luminosity dependences in B eq in our own sample, and 
in a combined sample where our sources are put together 
with radio-loud quasars from other studies. We discuss our 
findings in Section [], and the results are summarized in Sec- 
tion y. 

Unless otherwise specified, we have assumed an 
Einstein-deSitter universe model with Hq — 50 km s _1 
Mpc -1 and qo = 0.5 throughout the analysis. Also, in this 
paper, 'quasar' should be interpreted as 'radio-loud quasar'. 



2 QUANTIFYING GALAXY EXCESS 

2.1 Galaxy— quasar cross-correlation amplitude 

The use of the amplitude of the spatial cross-correlation 
function as a m easure of galaxy excess was initi ated b y Seld- 
ner & Peebles ( 197^) and Longair & Seldner ( |l97S| ). Since 
then, it has been widely used by other investigators for quan- 
tifying galaxy environments around various types of AGN 
at both lower and higher redshifts (e.g. Yee & Green 1984 
YG87; Prestage & Peacock 1988, 1989; Yates et al. 1989 
EYG91; Smith, O'Dea & Baum 1995; Wurtz et al. 1997 
De Robertis, Yee & Hayhoe 1998). We only give the expres- 
sions here, but a detailed derivation may be found in Longair 



& Seldner ( [19791 ). 

The spatial cross-correlation function, £ (r), is found to 
be well fit by a power-law with slope 7 and amplitude B eq , 

£(r) = J B gq r^, 

where B m is the amplitude we seek in this analysis. In order 
to determine _B gq from the data, the amplitude, A gq , of the 
angular cross-correlation function, u>(8) = AgqB 1 must 
first be determined. An estimator for this amplitude is: 

Atot - A b 3 - 7 T _i 



N h 



(1) 



which is directly obtainable from galaxy counts in the quasar 
fields and in images of background control fields. In the 
above expression, Mot is the number of galaxies in the 
quasar field within a circle of radius 9 centered on the quasar, 
corresponding to 0.5 Mpc at the quasar redshift. The num- 
ber of background galaxies from the control fields is denoted 
by Nh- The slope of the correlation function, 7, is usu ally 
assumed to be 1.77, as found by Groth & Peebles (1977) for 



the low-redshift galaxy-galaxy correlation function. YG87 
showed that the galaxy-quasar correlation fu nction has a 
slope consistent with this. Lilje & Efstathiou (1988), how- 
ever, found 7 = 2.2 in a study of the cross-correlation be- 
tween Abell clusters and the Lick galaxy counts down to 
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scales of 0.2 Mpc, but Prestage & Peacock (1988, 1989) have 
demonstrated that _B gq is insensitive to an incorrect choice 
of 7, provided that 7 ~ 2. We have therefore made the usual 
assumption that 7 = 1.77. 

The spatial correlation amplitude is obtained from 
A gq by the relation: 

AfgAgq 



B 



d. 



7-3 



where A/g is the average surface density of background galax- 
ies, and dg is the angular diameter distance to the quasar. 
The two quantities in the denominator are / 7 = 3.78, which 
is an integration constant, and $ (mn m , z), which denotes a 
universal luminosity function integrated down to the appar- 
ent completeness limit at the quasar redshift: $ (mn m , z) = 
ii(m z) ^ F ^^ le m t e S ra ted luminosity function gives 

the number of field galaxies per unit volume at the quasar 
redshift with magnitudes brighter than the completeness 
limit of the data. By normalizing B gq with $(mn m ,z), an 
absolute measure of the galaxy excess at the quasar red- 
shift is obtained, but unfortunately the luminosity function 
is poorly constrained at higher redshifts. This a critical step 
in the analysis, especially since B gq is inversely proportional 
to the normalization, (j>* , of this luminosit y func tion. As ad- 
vocated by YG87 and Yee & Lopez-Cruz ( |i9og| ), it is there- 
fore important that <f>* is consistent with the background 
galaxy counts obtained from the data. In our analysis, we 
have used the Can ada-France Red shift Survey (CFRS) lu- 
minosity function (Lilly et al. 1995) to track the evolution of 
the characteristic luminosity L* , and then used the observed 
galaxy counts in order to fit cj>* as a function of redshift. This 
is described in more detail in Section [| 

2.2 A more direct measurement 



The quantity A0.5, first introduced by Hill & Lilly ( f!99l| ), 
provides a more direct, but less deep measure of the cluster- 
ing strength than B eq . The galaxy counting is restricted to a 
certain magnitude interval, but a separate calculation of the 
luminosity function is not required. Hill & Lilly used it for 
quantifying galaxy environments around a sample of radio 
galaxies. For a radio galaxy with apparent magnitude m g , 
A?o.5 is the number of excess galaxies with magnitudes in the 
range m g to m g + 3 within a radius of 0.5 Mpc centred on 
the radio galaxy. Since the quasar appears as a bright stellar 
object, the magnitude of the host galaxy cannot be directly 
determined. Instead, we have estimated m g for the quasar 
host fro m the m agnitude-redshift relation for radio galaxies 
of Eales jl985|): 



(R) = 21.05 + 5.3 log; 



(2) 



valid for 19.0 < m g (R) < 21.3. Some of our quasar fields 
were also imaged in /, and in this case we estimated m g (I) 
by subtracting the R — I colour of an E/SO galaxy model 
(Guiderdoni & Rocca-Volmerange 1988) from m g (R). The 



iVn.5 measurement may suffer from systematic errors with 
redshift if radio galaxies and clusters galaxies evolve differ- 
ently in luminosity. At these redshifts, however, the effect is 
likely to be small, and we have therefore not taken it into 
account. 
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Figure 1. Radio luminosity at 408 MHz of quasars in the sample 
as a function of redshift. The luminosities were calculated from 
flux densities and spectral indices given in Table 1. 



3 SAMPLE SELECTION 

The sample investigated here consists of 21 radio-loud steep- 
spectrum (a > 0.5 where radio flux density is S v oc v~ a ) 
quasars with redshifts 0.5 < z < 0.82 covering the radio 
luminosity range 23.8 < log (L408MHZ /WHz^sr" 1 ) < 26.7. 
The sources were divided into high- and low radio lumi- 
nosity sub-samples selected from two different radio/optical 
flux-limited samples in order to better distinguish between 
luminosity and redshift dependences. 

The eleven sources in the low-luminosity sample were 
drawn from the 7C quasar (7CQ) survey (Riley et al. 1999) 
which consists of 7C sources with a flux density limit of 
Sisimhz > 0.1 Jy (« 0.05 Jy at 408 MHz) identified with 
objects on POSS-I (Palomar Observatory Sky Survey) E- 
plates with E < 20 and colour O - E < 1.8 (E ~ R 
and O ~ B). Consequently, the 7CQ survey finds only the 
optically brigh test s teep-spectrum quasars at any redshift. 
Willott et al. (1998) estimate the 7CQ survey to be 40 per 
cent complete due to the bright optical magnitude limit. 
However, at z = 0.5 and z = 0.8, the optical flux limit 
corresponds to Mb = —21.7 and Mb = —22.7 respectively 
(for AT-corrections, see below), hence according to the defi- 
nition employed by Willott et al., where a quasar must have 
Mb < —23, the 7CQ sample is complete in our redshift 
range. Clearly though, this is a rather artificial distinction, 
and we will be missing many broad-line objects with fainter 
nuclei. 

The eight sources in the high-luminosity sample were 
drawn from the Molonglo/APM Quasar Survey (MAQS) 
(Serjeant 1996; Maddox et al., in preparation; Serjeant et 
al., in preparation) which consists of sources with radio flux 
densities S408MHZ > 0.95 Jy from t he Molonglo Reference 
Catalogue (MRC) ( Large et al. 1981 ) identified with objects 
on the UK Schmidt plates having bj < 22.5. The MAQS 
thus have a faint optical completeness limit and a relatively 
bright radio flux density limit, and its comp letene ss was es- 
timated to be 99 per cent by Willott et al. jl998[ ). 

We also included two more quasars in the sample, 3C380 
and 5C6.189. The most powerful source, 3C380, was selected 
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Table 1. Basic data of the sources in the radio-loud quasar sample. The spectral index, a, is defi ned in the 
sense that r adio fl ux density is S y oc v~ a . The data on the MRC and 7C quasars are fro m Ser jeant (1996) and 



Riley et al. (1999), respectiv ely. R edshift and flux density of 5C6.189 are from Rossitter (1987) and the redshift 
of 3C380 from Laing et al. (1983). The flux densities at 408 MHz of the 7C quasars were calculated from flux 
densities at 151 MHz, assuming spectral indices as given in the table. 
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a The radio sizes and morphologies were determined from maps by Kapahi et al. (1998) (MRC 0032—203 and 
0106-233), Morganti, Killeen & Tadhunter (1993) CM RC 01 59-117 and 0405-123), Rossitter (1987) (5C6.189), 
Downes et al. (1986) (MRC 0222-008), Reid et al. ( |l995| ) (3C380) and Riley et al. (1999) (7C sources). An 
L-band VLA map of MRC 0406-180 is presented in Best, Rottgering & Lehnert (1999) where it is consistent 
with being unresolved, and it is used as a VLA phase calibrator at L-band in A-array, so must have an angular 
size <^0.5 arcsec. Radio sizes for MRC 0033—000 and 0144—058 are from Blundell (private communication). 
b 3C380 appears compact, but its intrinsic size is prob ably larger than 30 kpc since there is evidence for 



superluminal motion in this source (Wilkinson et al. 1990) 



from the 3CRR sample of Laing, Riley & Longair (1983), 
which has a flux limit of Si78MHz > 10.9 Jy (~ 6.05 Jy 
at 408 MHz). The radio quasar 5C6.189 was selected from 
t he 0. 1 < 5i5iMHz < 0.5 Jy complete sample of Rossitter 
(1987), which has essentially the same selection criteria as 



the 7C sources, so 5C6.189 is included in the low-luminosity 
sample. This quasar has a spectral index of 0.36, so strictly 
speaking it does not belong with the other steep-spectrum 
sources in the sample. It is unresolved to a limit of 3 arcsec 
in radio images, and since it also has a low radio luminosity 
we believe that it might be the Doppler boosted core of an 
intrinsically much less luminous radio source. We therefore 
neglected this object in the correlation analysis in Section Q 

The majority of the sources in the sample exhibit FRII 
morphology with bright cores, jets, edge-brightened lobes 
and hotspots, typical of radio-loud quasars. Two of the 7C 
and four of the MRC quasars are compact steep-spectrum 
(CSS) sources. CSS sources have apparent sub-galactic sizes 
(< 30 kpc) in radio and are thought to be intrinsically com- 
pact rather th an being extende d sources shortened by pro- 
jection effects ( Fanti et al. 199C ). It is believed that the CSS 
sources are either young radio sources in the stage of forma- 
tion, or that they are confined by a high-density ISM in the 
host galaxy which is likely to be denser than the IGM/ICM 



confining the more extended radio sources. It is uncertain 
whether this will bias the CSS sources in the sample toward 
richer or poorer galaxy environments, since it is difficult to 
see a reason for a correlation between the local gas density 
in the core of the host galaxy and the gas density (or the 
galaxy density) on Mpc scale in the surrounding cluster. 

In Fig. [l] we have plotted the radio luminosity of the 
quasars in the sample as a function of redshift, and various 
properties of the sources are listed in Table |l|. The quasar 
B luminosities in this table were calculated as follows: For 
the MRC quasars, we obtained tub from APM magnitudes 
assuming ttiapm ~ 29.42 — bj (Maddox, personal c ommuni- 
cation) and bj w m B ~ 0.09 (Metcalfe et al. 199 1|). For the 



7C sources, we calculated uib from E magnitudes by first 
assuming R~ E and thereafter converting to B by assum- 



ing a power law spectrum 



and using the B and 



R zero points given by Longair ( 1981 ) . The tub magnitudes 
of 5C6.189 and 3C3 80 we re found from R and B ma gnitudes 
given by Rossitter (1987) and Netzer et al. (1991;), respec- 
tively. The apparent magnitudes were thereafter converted 
to absolute Mb luminosities by assuming a if -correction 
equal to 2.5 (a — 1) log (1 + z). 

All sources, except 3C380 and 5C6.189, have galactic 
latitudes |6| > 42° and galactic reddening E (B - V) < 0.04. 
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Table 2. Positions and details about the observations of the quasar (upper part) and background fields (lower part). 
'NOT' is the Nordic Optical Telescope, 'McD' is the 107" telescope at the McDonald Observatory and 'HST' is the 
Hubble Space Telescope. 

Source RA (1950) DEC (1950) Filter Exp. time Telesc., date 
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Short exposure to enable deeper non-photometric images of this field to be properly calibrated. 



3C380 and 5C6.189 lie at latitudes 23?5 and 27?3 and have 
galactic reddening of 0.07 and 0.06, respectively. 



4 OBSERVATIONS 

Most of the quasar fields were imaged with the HiRAC 
camera at the 2.56-m Nordic Optical Telescope (NOT) at 
La Palma, Spain, during 1994 Dec 23-25 and 1997 May 10- 
15. The 7C3222 and 7C3450 fields were imaged with the 
107" telescope at the McDonald Observatory in Texas dur- 
ing 1996 Feb 17-22. Images were also taken of several back- 
ground control fields to provide a measurement of the field 
galaxy counts. The least biased way of obtaining background 
galaxy counts is to estimate it from the edges of the quasar 
fields, but the HiRAC CCD at the NOT has too small a field 
of view (3x3 arcmin) for this approach. The background im- 



ages were therefore obtained by offsetting the telescope 5 or 
10 arcmin either north or south of sources in our AGN sam- 
ple, or by selecting a random position in the sky at a similar 
galactic latitude as the programme targets. 

To extend our low radio luminosity sample, we used 
data on five 7C quasars obtained with the Hubble Space 
Telescope (Serjeant, Rawlings & Lacy 1997; Serjeant et 
al. 1999, in preparation). These images were taken in the 
F675W-filter in 4x500 s exposures with the WFPC2. In Ta- 
ble El positions and details about the observations of the 
quasar and background fields are provided. 

Eight of the quasar fields were imaged in two filters, 
either V and R, or R and /, depending on the redshift of the 
quasar, in order to straddle the 4000 A break in the quasar 
rest frame. This was done with the aim of using colours to 
study the evolution of the cluster galaxies, and to isolate 
cluster members with strong 4000 A breaks. These studies 
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Table 3. The observing runs. 



Date 


Tclcsc. 


CCD 


FOV 


Pixel scale 


Average 


Extinction 


Remarks 








(arcmin 2 ) 


(arcsec/pix) 


seeing (arcsec) 


in V (mag) 




1994 Dec 23,24 a 


NOT 


lk SiTe 


3x3 


0.176 


0.6 


0.13 


photometric 


1994 Dec 25 


NOT 


lk SiTe 


3x3 


0.176 


1.0 


0.13 


variable transparency 


1996 Feb 17-22 


McD 


lk TEK 


11x11 


0.700 


2-3 


0.17 


partly cloudy 


1996 July 21-23 


NOT 


lk SiTe 


3x3 


0.176 


0.7 


0.1-0.2 


photometric 


1996 July 24 


NOT 


lk SiTe 


3x3 


0.176 


0.7 


0.34 


dust in the air 


1997 May 10 6 -14 


NOT 


2k Loral 


3.7x3.7 


0.110 


0.7-0.8 


0.13 


some nights clear, 



some with cirrus 



a Most of the radio-loud quasar fields were obtained during Dec 23—24. 

b On May 10 the CCD had lost its quantum efficiency in blue resulting in a 0.5 mag reduction in the /-band sensitivity. 



will be discussed in a future paper. Quasars with redshifts 
z < 0.67 were imaged in V and R, whereas quasars with 
z > 0.67 were imaged in R and I. For the fields that were 
imaged in only one filter, R was chosen for z < 0.67 and I for 
z > 0.67. Six of the background fields were also imaged in 
two filters. With a few exceptions, the exposure time in every 
quasar and background field is 4x600 s (or 9x300 s), so that 
all fields have a similar depth. The images from the 107" 
telescope have exposure times of 11x300 s since the CCD 
was less sensitive than the CCD at the NOT, and the seeing 
was poorer. The weather and seeing conditions during the 
observations, as well as details about the CCDs that were 
used, are given in Table ^[ For standard star calibration, 
we observed fields from Christian et al. (1985) during all 
observing runs. 



5 DATA REDUCTION AND PHOTOMETRY 

The preprocessing of the images was performed with iraf. 
Twilight flats were used to flatfield the NOT data, but the 
/-band data from the 1997 May run were also flat-fielded 
with a skyflat formed by medianing 9-10 frames in order 
to remove fringing. After co-aligning the images, they were 
averaged together and cosmic rays were removed. 

The McDonald data were corrected for dark current and 
flat-fielded with either twilight or dome fiats. Since these 
data suffer from a partially cloudy sky, some extra reduction 
steps were found necessary. All images were divided with a 
skyflat in order to remove flat field variations during the 
night, and to improve the removal of the vignetting at the 
edge of the field. To take out the effect of clouds, some of 
the images were scaled so that some reference star in each 
image had the same flux as measured in the best image. 
Each scaled frame was then multiplied with a weight-factor 
inversely proportional to the variance before co-aligning and 
combining. 

The HST data were corrected for bad pixels using 
the bad pixel maps, and thereafter combined using the 



iraf/S"SDAS task CRREJ in order to remove cosmic rays. 



5.1 Photometry and object detection 

For object detection and photometry of faint galaxies in the 
CCD images we utilized FOCAS (Faint Object Classification 
and Analysis System, e.g. Valdes 1989). The detection limit 



in FOCAS was set to 2.5a above sky level and the 'built- 
in' detection filter was used, a 5x5 pixel array designed to 
increase the sensitivity to extended objects with low surface 
brightness. Furthermore, objects were detected only if they 
had sizes larger than a minimum size which was determined 
for each frame by measuring the average seeing from stars 
in the image. Objects were therefore recognized only if they 
had N contiguous pixels above the sky level, where N was 
varied to match the seeing and the pixel scale of the CCD. 
To estimate the magnitudes of the detected objects we used 
FOCAS total magnitudes. The total magnitudes in FOCAS are 
evaluated by filling in concavities in the isophote shapes until 
the area is twice the isophotal area and then measuring the 
flux above sky level within this expanded region. The FOCAS 
total magnitudes have been shown to give good estimates of 
true magnitudes (Koo, Ellis & Windhorst 1989). 



No attempt was made to separate stars and galaxies in 
the images, because the wide range in seeing conditions in 
which the data were taken made it hard to perform such 
separation reliably. Since the background fields were taken 
at similar galactic latitudes as the quasar fields, the num- 
ber of stars is approximately the same in both. Stars will 
therefore statistically be accounted for in the analysis. Also, 
at high galactic latitudes, galaxies are expected to dominate 
the counts fainter than R ~ 20. Each detected object was 
classified by FOCAS by using a template point spread func- 
tion, which was formed from typically 4-5 selected stars in 
each image. When we counted the number of objects in the 
images, we excluded objects that were classified as 'n' (noise) 
or 'long', and also in a few cases, objects coinciding with the 
spikes of bright stars. 

Magnitudes measured in the F675W-filter were cali- 
brated to Johnson R by using photometry of objects in fields 
that were imaged both with the HST and the NOT. With 
this approach, we found a zero point of 21.78±0.47 mag, 
which is similar to the one obtained by using the photflam 
and photzpt header keywords to compute a magnitude in 
the STMAG system, and then convert to Joh nson R as de- 



scribed in the WFPC2 Instrument Handbook (Biretta et al 



1996) 



All magnitudes were corrected for g alactic extinction 
using the maps of Burstein & Heiles (1982) and the G alactic 
extinction law by Cardelli, Clayton & Mathis ( 1989] ). 
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5.2 Completeness simulations 

To find the completeness limit and to gain as much informa- 
tion at faint magnitudes as possible, we simulated the loss 
of galaxies (in the ground-based data) in both the quasar 
and the background images due to incompleteness. We uti- 
lized ARTDATA in I RAF to add a grid of ~70 artificial galaxies 
with known magnitudes (typically 20-25 mag) and half-light 
radii distribution according to the HST Medium Deep Sur- 
vey (Schmidtke et al. 1997). Two different galaxy profiles 



been corrected for incompleteness as described in section 



were chosen in order to investigate which galaxies were most 
easily picked up by FOCAS, exponential discs resembling spi- 
ral galaxies and de Vaucouleurs profiles resembling elliptical 
galaxies. Inclination and position angles were randomly dis- 
tributed. 

After adding the artificial galaxies, the CCD frames 
were processed in FOCAS and the result of the detection 
process examined to find how many of the added galaxies 
were detected, and what their total magnitudes had been 
estimated to. This process was repeated several times, each 
time turning the flux of the galaxies down by 0.1 mag, until 
the fraction picked up by FOCAS had decreased significantly. 

By comparing the FOCAS-computed total magnitudes 
with the magnitudes originally assigned to the galaxies, we 
were able to estimate the errors made by FOCAS as a function 
of magnitude. These simulations show that the NOT data 
are complete down to 23.5 in R and 23.0 mag in / with 
errors of ±0.3 mag at the faintest limit. The data from the 
107" telescope are less deep with a completeness limit of 
22.5 in R and 21.5 mag in / with errors of ±0.4 mag at the 
limit. Furthermore, galaxies with exponential disc profiles 
were more easily detected than those with de Vaucouleurs 
profiles. 

In intervals of half a magnitude, we calculated the frac- 
tion of galaxies detected in FOCAS by averaging the detected 
fractions of exponential and de Vaucouleurs profiles. In or- 
der to correct for incompleteness at faint magnitudes, the 
number of galaxies in every magnitude bin was divided by 
the detected fraction in the corresponding bin, but only as 
long as the correction was less than a factor of 2. The cor- 
rected counts were used to find deeper and incompleteness- 
corrected values of the correlation amplitudes. The WFPC2 
images are deeper than the ground-based images in the sense 
that fainter galaxies will be more easily detected. However, 
since we only have ground-based background frames to nor- 
malize the counts, we took the same magnitude cut-off in 
the HST images as in the NOT images when computing 



6 RESULTS 

6.1 Field galaxy counts 

The background control images obtained at the NOT yield 
an area of 18 arcmin 2 in V, 58.7 arcmin 2 in R and 73.7 
arcmin 2 in 7, most of them taken during the observations in 
July 1996 and May 1997. The background images from the 
107" telescope cover an area of 103.3 and 97.7 arcmin 2 in R 
and I, respectively. Fig. |^ and [] show the number counts per 
deg 2 of field galaxies as obtained from the NOT images. The 
average number counts of galaxies in the quasar fields are 
also plotted. The number counts in these two figures have 



5.2, and the corrected field counts have slopes of 0.29±0.05 
and 0.34±0.06 in R and I, respectively. In the i?-band, the 
slope of the raw number counts is 0.24±0.03, whereas for the 
raw / counts, the slope is 0.26±0.06. Our estimates of the 
slopes are in agreement with that found by others (Smail et 
al. 1995; Steidel & Hamilton 1993; Lilly, Cowie & Gardner 
1991). For example, Smail et al. find a slope of 0.321±0.001 
in R and 0.271±0.009 in I, which is within the errors of the 
slopes of our completeness corrected counts. 

The error bars in Fig. ^ and [] were calculated as Poisson 
errors multiplied with a factor («1.3) that approximates the 
departure from Poisson stati stics if there is a non-ra ndom 
distribution of field galaxies ( Yee fc Lopez-Cruz 1999 ). 



6.2 

In the quasar fields, we counted galaxies with magnitudes 
brighter than the completeness limit within a circle of 0.5 
Mpc centered on the quasar (0.5 Mpc at z = 0.7 corresponds 
to ~ 66 arcsec). The average background galaxy count down 
to the same magnitude was thereafter subtracted in order to 
give the excess number of galaxies, N to t — Nb- Typical values 
of TVtot in the richest fields were 70-80 whereas the poorest 
fields had iVtot ~ 20-30. The background count, iVb, was 
typically around 45. A corrected excess was also calculated 
using the incompleteness corrected counts. In most fields 
this could be done down to a fainter magnitude than what 
we defined to be the completeness limit, but not fainter than 
24.5 in R and 23.5 mag in I, since here the correction factor 
was greater than 2. 

In order to prevent systematic errors, the quasar and 
background frames were treated identically during data re- 
duction and object detection. Also, the background galaxy 
counts obtained at one of the two (NOT and 107" ) different 
telescopes were applied only to quasar frames from the same 
telescope. The galaxy excess in the HST images was found 
by subtracting the NOT background counts. For fields that 
were imaged in two filters, we used the galaxy counts in the 
image taken through the reddest filter for computing -B gq , 
since this gives preference to red galaxies that are likely to 
lie at the quasar redshift. 

The field size of 3x3 arcmin at the NOT was sufficient 
to cover the 0.5 Mpc radius, provided the quasar was well 
centered on the frame. Two of the quasar frames were nev- 
ertheless found to have parts of the circle lying outside the 
edges of the field. No large areas were missing, 0.28 arcmin 2 
at the most, which is less than 8 per cent of the area of the 
circle. We also took out m 6.5 per cent of the area in the 
MRC0222-008 field because of a large, bright foreground 
galaxy. The missed and the chopped-out areas were cor- 
rected for by assuming a density of galaxies equal to the 
mean density within the remaining area. The same approach 
was applied to the HST data, where the 'L'-shaped field of 
view of the WFPC2 covered approximately half the area 
that we were interested in. 

The excess number of galaxies in each field was used 
to find the angular cross-correlation amplitude, A S q. Com- 
puted values of A gq , both from the raw counts and from the 
corrected counts, are shown in Table ||. We also averaged to- 
gether the galaxy number counts within the 0.5 Mpc radius 
in all the quasar fields and plotted it together with the field 
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Figure 2. This plot compares the average number count of galax- 
ies in the quasar (open circles) and the background images (filled 
circles). A clear excess of galaxies in the quasar fields is seen from 
R PS 20. The completeness limit is indicated by the vertical dotted 
line. The slope of the background counts at Ft > 20 (where the 
counts are not contaminated by stars) was found to be 0.29±0.05. 



galaxy counts as shown in Fig. ^ and |§] to demonstrate that 
there is a real excess of galaxies. 

In order to obtain the spatial cross-correlation ampli- 
tude, -Bgq, the angular correlation amplitude has to be nor- 
malized with the number of field galaxies at the quasar 
redshift, as given by the integrated lumi nosity function , 
& (miim, z ). We used the Schechter function ( [Schechter 197!: ) 
with parameters a — — 0.89 and M* AB (B) — —20.83 based 



on Lilly et al.'s ( 1995 ) estimation of the CFRS luminos- 
ity function at redshifts 0.5 < z < 0.75. The values of <f>* 
employed were 0.0072 and 0.0052 Mpc -3 in R and I, re- 
spectively. Below, we explain how the <f>*'s were chosen in 
order to be consistent with the observed field galaxy counts. 
The CFRS luminosity function is well suited to our analysis, 
especially since CFR S is a n /-band selected survey, and re- 
sults from Lin et al. (1999) based on the CNOC2 (Canadian 
Network for Observational Cosmology) Field Galaxy Red- 
shift S urvey show that there is general agreement between 



the CF IS and the CNOC2 luminosity function. 



To make sure that the observed background galaxy 



counts were consistent with the </>* normalization, we con- 



structed galaxy counts from a combination of four Schechter 
function segments. These constructed counts were compared 
to the observed counts to find the <f>* that best matched the 
data. We selected luminosity function segments for four red- 
shift intervals, (0.0, 0.2), (0.2, 0.5), (0.5, 0.75) and (0.75, 1.0). 
In the first interva l, we u sed the Schechter parameters found 
by Loveday et al. ( 199S| ) for the local luminosity function in 
the Stromlo-APM survey; M bj = -19.5, a = -0.97 and 
(f>* = 0.014 Mpc -3 . For the three other intervals, we used 
M% B (B) = -21.04, -20.83 and -21.24 as found by Lilly 
et al. for their 'AH' colour sample in the CFRS. We also 
fixed a at —0.89 in these three redshift intervals, since the 
slopes reported by Lilly et al. at these redshifts are seen to 
vary a lot, which presumably is not a real evolutionary ef- 
fect. Since the CFRS luminosity function is defined over a 
narrow magnitude range, the faint-end slope may have high 



Figure 3. Same as in Fig. 2, but for counts in the /-band. The 
slope of the background counts at / > 20 is 0.34±0.06. The 
'bump' at / < 20 is probably caused by stars. 



uncertainties. For example, in the redshift range (0.5, 0.75) 
that we are most interested in here, and at luminosities from 
Mab (B) = -19.7 to -22.9, Lilly et al. find that a = -0.5. 
The faintest absolute magnitude in this range corresponds 
to w —21 in R, whereas our completeness limit in R corre- 
sponds to « —20 at z — 0.6, so our data are approximately 
one magnitude deeper. We therefore chose a = —0.89 since 
this is the slope they find for the total luminosity function 
from z — 0.0 to 1.3. 

The four luminosity function segments were integrated 
with respect to luminosity down to the completeness limit 
of our data at the average redshift in each interval. The 
completeness magnitudes in R and I were found by applying 
K- correct ions and galaxy colours of an Sa galaxy listed as 
a, func tion of redshift by Rocca-Volmerange & Guiderdoni 
and Guiderdoni & Rocca-Volmerange (1988). The 



characteristic magnitudes in each interval were transformed 
f rom Bab to R and / by first using that Bab = B — 0.17 
( |Oke 1972[ ), and thereafter assuming B — V = 0.74, V - 
R = 0.68 and R — I — 0. 57 which are the colours of an Sa 
galax y model at z = (Guiderdoni & Rocca-Volmerange 



198£). The M bj magnitude in the lowest reds hift interval 
was conv erted to B by assuming bj — B - 

1)- 



al. 1991 



0.09 (Metcalfe et 



Then the integration was carried out over redshift, from 
z = 0.0 to 1.0. From these approximate calculations we con- 
structed number counts that were subsequently matched to 
the observed counts by varying the 0*'s of each segment 
until a minimum in \ 2 was reached. For the /-band data 
we found it sufficient to use the same scaling for all inter- 
vals, and the best fit was obtained with (j>* = 0.0052 Mpc -3 
for the (0.5, 0.75) interval. In R, very little adjustment was 
needed in 0* to fit the data, and for the (0.5, 0.75) interval 
we found a best fit with 0* = 0.0072 Mpc" 3 . 

To normalize B sq , we evaluated $(mu m , z) for every 
quasar field by integrating the Schechter function from — oo 
to an absolute magnitude corresponding to the complete- 
ness limit at the quasar redshift. The absolute limiting mag- 
nitude in every quasar frame was found by converting the 
apparent completeness limit to an absolute magnitude using 
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Jf-corrections for a hot E/SO galaxy ( Rocca-Volmerange &j 



Guiden loni 



Typical values obtained for the integrated 



luminosity function were $ (roiin 



0.005-0.007 Mpc" 



The results of the computations of _B gq are shown in 
Table |^. The errors in A eq and B gq were calculated according 
to 



[(iVtot - JV b ) + 1.3 2 7Vt 



I 1/2 



ATt, 



(•3) 



( Yee & Lopez-Cruz 1999) and reflect the errors in the galaxy 
counts on the quasar and background images, but not un- 
certainties in the luminosity function. (Yee & Lopez-Cruz 
also argue that the inverse of this error underestimates the 
statistical significance of excess galaxy counts which may 
be better represented by (jV to t - N b )/1.3^N^). In Table ^ 
Bgq found from both the raw and the corrected counts are 
shown. On average the corrected B gq 's are larger since these 
were deduced from deeper counts, in some cases down to 
R = 24.5 (for MRC0144-058, MRC0159-117 and 7C3066). 
In cases where the quasar and background fields had dif- 
ferent completeness limits, we took the brightest magnitude 
as the limit. Some fields could not be pushed further with 
regard to the limiting magnitude, so these are left as 'un- 
changed'. The HST images were also not corrected. 

The counting of galaxies should not be done either too 
far below or too far above M* of the surrounding galaxy 
cluster in order to obtain a reliable result. At faint magni- 
tudes the background counts rise faster than the associated 
galaxy counts, and a possibl e exce ss may drown in the back- 
ground. Yee & Lopez-Cruz (1999) showed that if the appro- 
priate luminosity function is used to normalize the counts, 
i.e. that the <j>* of this luminosity function can reproduce 
the observed background counts, then the best results are 
obtained by counting down to a magnitude between M* + 1 
and M* + 3. As far as possible, we have attempted to reach 
M* + % and our completeness magnitudes nf R — 23 5 and 
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Figure 4. The cross-correlation amplitudes from this work, de- 
rived from the corrected counts, plotted together with ZJ gq for 
steep spectrum radio-loud quasars from the literature. The 
galaxy-galaxy correlation amplitude for low-2 galaxies in the 
field, S gg fa 60 Mpc 1 - 77 (Davis & Peebles 1983), is shown as 
a dotted line across the plot. 



to M* if the quasar fields are contaminated by foreground 
galaxies. This approach also assumes that the field and clus- 
ter luminosity function evolve together. Second, they scaled 
(f>* of the evolved local luminosity function to match the 
observed background counts. 

For (z) = 0.6, YG87 used M* = -21.55, a = -1.2 and 
<t>* = 0.0058 Mpc" 3 (50 = 0.5), and EYG91 used M* = 
-21.73, a = -1.0 and 0* = 0.0027 Mpc" 3 (go = 0.02). 
To compare the characteristic magnitudes, we transformed 
from r to R by assumin g r — R ~ 0.43 (Fukugita, Shi- 



/ = 23l0 correspond to Mb 



-20.5 and M, 



-20.8 at 



masaku fc Ichikawa 1995) and applied -corrections for a 



hot E/S0 galaxy ( jGuiderdoni fc Rocca-Volmerange 198£ 
since YG87 and EYG91 give M* in the observed wave- 
band. We found M R (YG87) = -22.98 and M* R (EYG91) = 
-22.85 (g = 0.5). By adding a colour of R - I = 0.78 to 
the R magnitudes, we obtained M| (YG87) = -23.76 and 
Mj (EYG91) = -23.63. It therefore seems that we have 
used an M* that is approximately one mag fainter than in 
both of these studies, and also a <f>* that is larger. A fainter 
characteristic magnitude will lower the value of the inte- 
grated luminosity function, thereby bringing B gq up by a 
small amount, but at the same time, a larger value of <j>* 
will increase the number of galaxies above the limiting mag- 
nitude, with the result that B gq decreases. 

The mean B gq for the quasar fields in our sample is 
265 Mpc 1,77 with an error in the mean of 74 Mpc 1,77 (from 
a combination of intrinsic dispersion in £? gq and measure- 
ment error). With 5C6.189 excluded, this becomes 249T75 
Mpc 1 ' 77 . To see how the choice of different luminosity func- 
tions affected our measurements, we recomputed B gq with 
both YG87's and EYG91's luminosity functions. In the first 
case, we obtained overall lower values of B sq , the mean was 
found to be 185±57 Mpc 1,77 . In the second case, we trans- 
formed to a go = 0.02 cosmology for consistency, and found 
a better agreement, the mean B gq was 257±77 Mpc 177 . The 
correlation amplitude is rather insensitive to go, and using 



redshifts 0.6 and 0.7 (assuming A"-corrections of approxi- 
mately 1 and 0.5 mag). Provided that M* of the associated 
galaxies is approxim ately the same as the one in the field 
given by Lilly et al. (1995), we have reached deep enough in 



all fields, except 7C3222 and 7C3450. These two fields have 
shallower completeness limits (imaged at the 107" telescope) 
of 22.5 in R and 21.5 in /, corresponding to M R m —21.5 
and Mj w -22.4, i.e. not as deep as M* + 1. Using the 
corrected counts, we go half a magnitude deeper and should 
therefore be on the faint side of M* + 1. 



6.3 Comparison with previous work 

We may compare with YG87 and EYG91 who also quanti- 
fied quasar environments in terms of B gq . The greatest dif- 
ference between their and our analysis is the method used 
for choosing the appropriate luminosity function for nor- 
malizing -B gq . Since determinations of the field galaxy lu- 
minosity function at higher redshifts were not available to 
them, they used local luminosity functions from the liter- 
ature which they fitted to the observations in two steps. 
First, they added an evolution in M* that matched the ob- 
served luminosity function of the associated galaxies. A pos- 
sible risk with this is that too much evolution may be added 
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Table 4. Clustering statistics for the quasar fields consisting of the angular and the spatial galaxy-quasar 
cross-correlation amplitudes, A gq and B gq , and the more direct Abell-type measurement, No. 5. The cor- 
relation amplitudes were calculated both from raw galaxy counts and from galaxy counts corrected for 
incompleteness. The magnitude limits were set to _R = 23.5 and I = 23.0 for the raw counts, whereas a 
deeper limit was applied for the corrected counts, maximum 24.5 in R and 23.5 in /. The fields that are 
denoted 'unchanged' are fields that could not be pushed above the completeness limits of 23.5 in R and 23.0 
mag in /. The corrected counts were not applied to the HST images either, as denoted by The errors in 
A gq and B gq are Ayl gq and A£? gq from Eq. 3. 
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cording to the definition of Kellerman et al. (1989). These 



smaller by a factor of ~ 1.15 if qo = 0.02 instead of 0.5. 

We have one field in common with Yee & Ellingson 



( 1993 ) with which we may compare more directly. They 
report that B gq in the MRC0405-123 field is 905+277, 
whereas we find 875+176 Mpc 1 ' 77 . When we recalculated 
B gq with EYG91's luminosity function, we found 1006+231 
for this field, and with YG87's luminosity function as nor- 
malization, we obtained 647+96 Mpc 177 . There thus seems 
to be very good agreement between our and the Yee & 
Ellingson result for this field, and in particular that this 
quasar lies in a relatively rich cluster. 

In Fig. ^ we have plotted our B gq measurements ob- 
tained from the corrected counts together with B gq for 
YG87's and EYG91's quasar fields as a function of redshift. 
In this figure, we have also added quasar fields reported by 
Yee & Ellingson ( 1993 ) that do not overlap with sources from 
YG87 and EYG91 (referenced in tab le 1 in Yee & Ellingson 



as coming from Yee & Green (1984) and from unpublished 
CTIO and CFHT images). We used photometric informa- 
tion available in the NASA/IPAC Extragalactic Database 
(NED) to find the spectral index of the literature sources, 
and found that some of them were flat-spectrum quasars. 
These were omitted in Fig. ^, since we only want to com- 
pare with other steep-spectrum sources. (We also found that 
three of the quasars, 0208-018, 0438-165 and 0449-183, 
listed as radio-loud in table 2 of EYG91, are radio-quiet ac- 



were also ignored.) It seems that there are no systematic dif- 
ferences in our estimates of £? gq and the other estimates from 
the literature. The majority of the literature quasar fields 
have amplitudes below 500 Mpc 1 ' 77 , a trend which is present 
also in our fields, and the overall impression from Fig. ^| is 
that the correlation amplitudes from the different samples 
are consistent with each other and fall approximately in the 
same range. We do not see any redshift dependence in -B gq , 
either in our own sample or in the combined sample, as 
shown in Fig. ^. Adding the flat-spectrum quasars seems to 
make no substantial difference to the plot, and it is difficult 
to see any significant trend in the flat-spectrum population 
alone. 



6.4 Comparison with other two— point correlation 
functions 

The average £? gq for our sample, 265+74 Mpc 1 ' 77 , is signifi- 
cantly larger than the galaxy-galaxy correlation amplitude 
in the field measured both locally and at higher redshifts. 
The local galaxy-galaxy correlation function, £ gg (r) = 
( r / r o) 1 , has been measured by several investigators (e.g. 
Davis & Peebles 1983; Loveday et al. 1995; Guzzo et al. 1997) 
and results emerging from these studies show that ro ~ 10 

= 50 



Mpc and 7 ~ 1.7, i.e. that B g 



60 Mpc 1 ' 77 (Hq 
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Figure 5. B gq plotted versus No. 5. The solid and the dotted 
lines show the best fit of B gq = (37.8 ± 10.9)^0.5- The dashed 
line is the relation for radio galaxies from Hill & Lilly (1991), 
B gg =30iVo.5. 
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Figur e 6. No. 5 for the sample of radio galaxies by Hill & Lilly 
(1991) (open circles) plotted together with -/V0.5 for the quasars 
in our sample (filled circles). 



km s _1 Mpc -1 ). Written as a correlation length, our mean 
Bgq corresponds to ro = 23.4 Mpc. 

Several attempts have been made to measure £ gg up to 
z ~ 0.5, but the resul ts so far seem to disagree to some ex- 
tent. Le Fevre et al. (1996), using the CFRS, find evidence 



for a decrease in the correlation length toward higher red- 
shifts wit h ro = 2.66 Mpc and 7 = 1.64 at z = 0.53. Hudon 
& Lilly (|1996|) find r = 3. 78 M pc at z = 0.48. O n th e 



other hand, Carlberg et al. (1998) and Small et al. (199£) 
find no evidence for evolution between z « 0.1 and 0.6, 
and Small et al. get a correlation length of ro = 7.4 Mpc 
at z ~ 0.3. Nevertheless, the values we are measuring for 
Bgq are larger than the observed galaxy-galaxy correlation 
amplitude at both higher and lower redshifts, and show that 
we have clearly detected enhancements in the galaxy density 
around the quasars. 

We can also compare our B g 



plitude of the clust er-ga laxy cross-correlation function, £ cg . 
Lilje & Efstathiou (1988) measured £ cg on scales down to 0.2 



Mpc using Abell clusters and the Lick galaxy counts. They 
found that £ cg was well fit by the same functional form as 
£ gg , but with 7 = 2.2 and a correlation length rp = 17.6 
Mpc. Subsequently, Croft, Dalton & Efstathiou ( 1999] ) have 



shown that this is also a good approximation to £ cg obtained 
from the independent Stromlo-APM galaxy and APM clus- 
ter redshift surveys on scales of 0.2 to > 20 Mpc. Prestage 
& Peacock (19 88, 19 89), Andersen & Owen ( |l994| ) and Yee 
& Lopez-Cruz ( 1999 ) have examined B cg measured in 1, 1- 
1.5 and 0.5 Mpc circles, respectively, as a function of cluster 
richness, all assuming 7 = 1.77. Comparison of these mea- 
surements of B cg with the one of Lilje & Efstathiou reveals 
some discrepancies, presumably caused by systematic dif- 
ferences in the measurements and scatter in the richnesses 
of the clusters observed. To compare the Lilje & Efstathiou 
function with the rest, we evaluated £ cg at a radius of 0.5 
Mpc, then converted to B cg assuming 7 = 1.77. This comes 
out to B cg — 740 Mpc 1 ' 77 . Prestage & Peacock estimate 
Beg « 690 Mpc 1,77 for an Abell class 1 cluster, and Ander- 
son & Owen « 615 Mpc 1,77 , both broadly consistent with 
the Lilje & Efstathiou function (which was calculated for 
clusters with Abell richness class > 1). In contrast, Yee 
& Lopez-Cruz using CCD data find a median B cg « 950 
Mpc 1 ' 77 for an Abell class 1 cluster. We have decided to 
take the mean for Abell class 1 to be 740 Mpc 1 ' 77 simply 
on the basis that more clusters were involved in the Lilje & 
Efstathiou study, and note that this is more likely to be an 
overestimate than an underestimate as clusters richer than 
Abell class 1 were included. We have two quasar fields that 
have Bgq larger than this value, which suggests that a few 
quasars at 0.5 < z < 0.8 are found in clusters of approxi- 
mately Abell class 1 or greater, but that most are in poorer 
groups or clusters. 

We believe that we have not overestimated the richness, 
since we also find fields with negative B gq . However, a small 
or negative B gq need not mean that there are no galax- 
ies associated with the quasar, but simply that the associ- 
ated galaxies are lost in the fluctuations of the background 
counts. 



6.5 Calculation of No. 5 



We also calculated N0.5 (Hill & Lilly 1991) for each quasar 
field by counting the number of galaxies within the 0.5 Mpc 
radius with magnitudes ranging from m g to m g + 3, where 
m g is the typical magnitude of a radio galaxy (Eq. ^) at 
the quasar redshift. The average number of galaxies within 
the same magnitude range in the background control im- 
ages was subtracted from the counts in the quasar fields. 
The estimates of No. 5 obtained in this manner are listed in 
Table § in the last column. Since the 7C3222 and 7C3450 
fields from the 107" telescope at the McDonald Observatory 
are not deep enough for m g + 3 to be brighter than the com- 
pleteness limit, we have calculated No. 5 from m g down to 
the completeness limit (21.5 in I and 22.5 in R). 

Although this is a measurement intended for fields 
around radio galaxies, it seems to give a good estimate of 
the richness in quasar fields since we find an overall cor- 
respondence between B gq and N0.5, as shown in Fig. |[ In 
this figure, large values of B gq is seen to pair rather well 
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with large numbers of No. 5, and fitting a straight line gives 
a best fit of B gq = (37.8 ± f0.9) N0.5, slightly steeper than 
the relation, B gg = 30iVo.5, found by Hill & Lilly for ra- 
dio galaxies, but still a good agreement considering that we 
have only estimated m g from a magnitude-redshift relation. 
There are a couple of fields with higher No. 5 that do not fol- 
low the trend as nicely as the rest of the fields (7C3814 and 
MRC 0144—058). Some fluctuations are expected, as N0.5 is 
not as deep a measurement as B gq , and also the number of 
field galaxies in the quasar images fluctuates to some level. 

In Fig. |(| we plot iVo.s for quasars in our sample to gethe r 
with N0.5 obtained for radio galaxies by Hill & Lilly (1991). 
(Hill & Lilly also assume Ho = 50 km s _1 Mpc -1 and go = 
0.5). Our values of N0.5 seem to lie approximately in the 
same range as those obtained by Hill & Lilly, although at 
the higher redshift end, Hill & Lilly's values for a few sources 
are larger. 



7 RADIO LUMINOSITY AND REDSHIFT 
DEPENDENCES OF THE CLUSTERING 
STRENGTH 

In order to examine whether the strength of galaxy clus- 
tering around the quasars is determined by epoch or radio 
luminosity, we looked for correlations between B gq , redshift 
and radio luminosity. 

In Fig. [j] we have plotted B gq as a function of radio 
luminosity at 408 MHz, I/4osmhz- The figure suggests that a 
positive correlation between B gq and radio luminosity may 
be present, at least there seems to be more of a trend with ra- 
dio luminosity than with redshift if we compare with Fig. ^ 

Our sample is rather small and covers a relatively nar- 
row redshift range, we therefore included the measurements 
of quasar environment s dow n to z as 0.2 from YG87, EYG91 
and Yee & Ellingson (1993). We used photometry from the 
NED to find flux densities and spectral indices to calcu- 
late radio luminosities for the literature sources. In a few 
cases where the source was not availa ble in NED, we use d 
the NRAO VLA Sky Survey (NVSS) ( |Condon et al. 1998| ). 
Since one of the aims was to investigate radio luminosity 
dependences, we excluded sources that were flat-spectrum 
quasars (nine sources in EYG91 and five from Yee & Elling- 
son, as well as 5C6.189 from our own sample), leaving us 
with a total of 71 radio-loud steep-spectrum quasars. (The 
three sources in table 2 of EYG91 found to be radio-quiet 
were also omitted.) 

In Fig. ^, we have plotted the B gq 's for all sources as a 
function of radio luminosity. To examine whether there are 
any trends in the data, we divided the quasars into redshift 



and luminosity sub-samples and computed the mean Bg q in 
each sub-sample. The result of this is shown in Table H. 

The mean B gq for all quasar fields in our sample, 
249±75, is seen to be consistent with the mean of 304±51 
for the literature sources, and the mean of 289±42 Mpc 1,77 
for the combined samples. The fact that the normalization of 
B gq was done in a self-consistent manner, both by us and the 
other investigators, is reassuring when we utilize their results 
in analyses like this. We also bear in mind the good agree- 
ment about the richness of the MRC0405— 123 field. As seen 
in Table [| the literature samples show a higher (B gq ) in the 
z > 0.5 sub-sample than in the lower redshift sub-sample, 
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Figure 7. B gq plotted as a function of radio luminosity at 408 
MHz. 



and we interpret this as the epoch-dependent evolution in 



B gq claimed by Yee & Ellingson (1993). The mean B gq for 



the combined samples is also seen to be larger in the high- 
redshift subsample, but here the difference smaller. However, 
the difference in the mean B gq for the low and high radio 
luminosity sub-samples is more pronounced. Here, we found 
(_B gq ) = 136 ± 72 at low luminosities and (B gq ) = 353 ± 49 
Mpc 1,77 at high luminosities. The sources in our own sam- 
ple gave a similar result, but with an even larger difference, 
(Bgq) = 107 ± 112 at low luminosities and 365±89 Mpc 1 ' 77 
at high luminosities. The division between the high- and low 
radio luminosity sub-samples was set at L408MHZ = 10 24 ' 75 
W Hz -1 sr _1 , which corresponds to the break in the radio 
luminosity function. 

In order to quantify the trends better, and to disentan- 
gle the effects that redshift and radio luminosity may have 
on B gq , we followed the approach of Yates et al. (1989) and 
computed the Spearman partial rank correlation coefficient, 
tax,y, for the different samples. The vax,y statistic is de- 
fined as 



tax,y 



TAX — TXYTYA 



[(1 



(1 



,1/2 > 



where tax , rxY and Tya are the usual Spearman rank corre- 
lation coefficients. The tax,y statistic is a number between 
— 1 and +1 and serves as an indicator of the trend of large 
values of A to be paired with large values of X when the 
effect of the third variable, Y, has been removed. A value of 
(— ) + l indicates a perfect (anti)correlation between A and 
X at constant Y. For N data points, the significance level 
of r A x,Y is 



Dax.y = -(iV-4) 1/2 ln 



1 + TAX,Y 



1 



- TAX,Y 



(Macklin 1982), and is approximately normally distributed 



about zero with unit variance if the null hypothesis, that 
the A-X correlation arises entirely from those of Y with 
A and X separately, is true. Here, we let A, X and Y be 
combinations of B gq , L408MHZ and z, so that the correlation 
between e.g. B gq and Z/408MHZ with z held constant can be 
found. By using this statistic, we hope to determine whether 
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Table 5. Mean B gq for sub-samples of steep-spectrum quasars. 
The results for our sample, the literature samples (YG87; EYG91; 
Yee &; Ellingson 1993) and all samples combined, are shown. Ra- 
dio luminosity, L408MHZ is given in W Hz -1 sr _1 and the unit 
of Bgq is Mpc 1 - 77 . The number of objects in each sub-sample is 
denoted by N. 
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Literature 



Our + 
literature 
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249±75 


20 


304±51 


51 


289±42 


71 


z <0.5 
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Table 6. Spearman partial rank correlation coefficients for the 
three different samples. Our sample, samples from the literature 
(YG87; EYG91; Yee & Ellingson 1993), and our sample com- 
bined with the literature samples. The lower half of the table 
shows the results of the correlation analysis when the compact 
steep-spectrum sources were excluded. The numbers in parenthe- 
sis behind the correlation coefficients are significances. 



Sample 



N 



IB 



gq^408> z 



r B„ a z,LA 



>-L, 



Our 20 0.41 (1.7<t) 

Literature 51 0.26 (1.8 a) 
Combined 71 0.39 (3.4a) 



0.02 (O.lcr) 0.22 (0.9cr) 
0.19 (1.3<t) 0.39 (2.8a) 
0.10 (0.8<t) 0.23 (1.9o-) 



Our 14 0.60 (2.2a-) 

Literature 44 0.27 (1.8<r) 
Combined 58 0.47 (3.7a) 



-0.18 (-0.6a) 0.41 (1.4<r) 
0.14 (0.9a) 0.39 (2.6a) 
0.02 (O.lcr) 0.22 (1.6a) 



there exists a correlation in _B gq with either z (independent 
of L408MHz) or L408MHZ (independent of z). 

The Spearman partial rank correlation coefficients were 
computed in three different cases: Our own sample, the lit- 
erature samples (YG87+EYG9f+Yee & Ellingson) and the 
samples combined. The resulting correlation coefficients and 
their significances are shown in Table |(| The table is divided 
horizontally, where the upper part shows the Spearman co- 
efficients calculated by including all steep-spectrum quasars, 
whereas in the lower part we have included only those steep- 
spectrum sources which have extended radio morphologies, 
i.e. we have excluded quasars that have sizes r < 30 kpc 
and are classified as CSS sources. The CSS sources may not 
sample the large-scale density distribution outside the ISM 
of the host galaxy, and therefore might not be expected to 
follow the same trends as the rest of the source population. 
The Spearman statistic without the CSS sources remains 
almost unchanged, however. 

The numbers in Table |(j seem to indicate a weak, but 
significant, dependence of £> gq on radio luminosity. Both in 
our sample and in the combined samples, it is the most sig- 
nificant correlation, and it is strongest when all the samples 



are combined, giving rs z, 4 



0.47 with a 3.7cr signifi- 



cance when the CSS sources are excluded. In the literature 
sample, it is the correlation between radio luminosity and 
redshift that dominates. This makes it difficult to interpret 
any changes in (B gq ) within the two sub-samples, and espe- 
cially we cannot be certain that there is a true evolution with 



Figure 8. B gq versus radio luminosity for objects in our sample 
together with steep-spectrum quasars from the literature. 



redshift. However, both in our sample and in the combined 
samples r_g gq z,L 40S is negligible and insignificant compared 
to fB gq L 4 o 8 ,z- The result of this analysis is therefore sugges- 
tive of a radio luminosity dependence in B gq , and we believe 
that there is no evidence for a redshift depende nce as has 
been claimed by e.g. YG87 and Yee & Ellingson (|l993j). 



We may consider that if our estimates of B gq are sys- 
tematically different from the others, then a bias could be 
introduced, but we do not think that this is the case. First, 
the agreement about the value of B gq for the field around 
MRC0405— 123 is very good. Second, the trends within our 
own sample are the same as those found for the combined 
samples. Third, using the B gq estimates determined from 
the raw counts in our data also give a significant correlation 
between B sq and L408MHZ, and produces the same trend as 
above when we combine with the literature sources. Fourth, 
the result of the correlation analysis using the recomputed 
B gq 's with EYG91's luminosity function as normalization 



However, we may note, as mentioned by Yates et al. ( 1989 ), 
that rB scl z,L 408 computed for a combination of samples that 
have been analysed by different authors should be viewed 
with some caution since differences in 3> (miim, z) and the 
assumed cosmology may introduce spurious effects. 



8 DISCUSSION 

8.1 Comparison with radio galaxy environments 

Previously, studies of radio galaxies spanning a wide range 
in radio luminosity at a given redshift have found that the 
richness of radio galaxy environments is increasing with red- 
shift, but does not correlate with radio luminosity (Hill & 
Lilly 1991; Arlington-Smith et al. 1993). It would be intrigu- 
ing if a correlation between clustering strength and radio lu- 
minosity was found for radio-loud quasars but not for radio 
galaxies. The 'Unified Scheme' for radio-loud AGN predicts 
that these are the same type of object viewed from different 
angles, so the richness of the environments should have the 
same behaviour with radio luminosity and redshift. Why, 
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then, are the radio quasars apparently different to the radio 
galaxies? 

First, we consider the possibility that selection effects 
are responsible. Our sample was initially chosen such as 
to consist of one low- and one high radio luminosity sub- 
sample. The sources in these two sub-samples were selected 
r andomly from th e complete radio/optic al flux-limited 7CQ 
flRiley et al. 1999[ ) and MAQS samples jSerjeant 1996| ), on 
the basis of redshift, right ascension and declination. The 
only possible source of bias would be a selection effect which 
leads us to give preference to less dense environments at low 
radio luminosities (the 7CQ quasars) and more dense envi- 
ronments at high radio luminosities (the MAQS quasars). 

This selection effect is possible if the radio luminosity 
of a quasar is controlled both by the environmental rich- 
ness/density (-Bgq) and the bulk kinetic power in the radio 
jets (op tica l luminosity). We discuss this in greater detail in 
section S.3. The 7CQ sample has a relatively bright optical 
flux limit, which selects sources with unusually high optical 
luminosities for their radio luminosities. The objects that we 
miss in the 7CQ survey are therefore optically faint with low 
jet powers, so their radio luminosity must be enhanced by a 
dense environment. We may therefore be biased against opti- 
cally faint quasars in dense environments and favour quasars 
in low-density environments instead, characterized by low 

In the MAQS sample, where the optical flux limit is 
faint, the effect would be to favour quasars in dense en- 
vironments, i.e. high -B gq 's. A faint quasar with weak jets 
in a dense environment can have the same radio luminos- 
ity as an optically bright quasar in less dense environments 
since the powerful jets of the latter will enhance the radio 
luminosity. Due to the steep increase of the optical quasar 
luminosity function towards low luminosities, we have more 
optically faint sources than we have bright sources, and since 
the sample also is radio flux-limited, we miss optically faint 
sources in low-density environments. The preference is thus 
given to optically faint sources in high-density environments 
having high B gq 's. 

In the YG87 and EYG9I samples similar effects may 
operate, as most radio quasar samples are selected with an 
optical flux limit close to the POSS-I plate limit, but they 
are harder to quantify as the samples are not complete. 

A second possible explanation is that we are not com- 
paring like with like, especially at radio luminosities close 
to the break in t he radio luminosity function. Recently, 
Laing et al. (1994) found evidence for two different types 
of radio galaxies: low-excitation narrow-line radio galaxies 
(LEG), appearing as typical radio galaxies in radio, but with 
weak optical emission lines that commonly have low excita- 
tion, and high-excitation narrow-line radio galaxies (HEG) 
with strong, high-ionization emission lines. The HEG's are 
thought to be the objects which unify with quasars, whereas 
the LEG's, unless seen with their jets pointing almost di- 
rectly along the line of sight (when they are seen as BLLac's) 
are radio galaxies at all orientations since no broad lines are 
observed. If the clustering properties of HEG's and LEG's 
differ, then a sample of radio galaxies with a mix of both 
HEG's and LEG's (such as that of Hill & Lilly 1991) must 
be considered less uniform and difficult to compare with 
quasars. 



8.2 The clustering properties of quasars and radio 
galaxies 

Because our quasars are typically in poor clusters, we can 
estimate their clustering (i.e. B qq ) by assuming the ratio 
B qq / B qc = B cc /B cg ~ 1.7, where B cc and B cg are the 
cluster-cluster and cluster-gal axy c orrelation amplitudes re- 
spectively, from Croft et al. (|l997j) and Croft et al. (1999). 



This gives a correlation amplitude for radio-loud quasars of 
~ 490 Mpc 1,77 (after correcting all numbers to a 7 = 1.77 at 
0.5 Mpc for consistency). This is large compared to lo w red - 
shift radio galaxies measured by Peacock & Nicholson ( 1991 ) 
who obtain -B g * g * ~ 240 Mpc 1,8 , and the faint radio g alaxie s 
in the FIRST survey, for which Magliocchetti et al. jl99^ ) 
estimate B g * g * « 220 Mpc 1 ' 8 . Both the Peacock & Nichol- 
son and the Magliocchetti et al. studies are dominated by 
lower luminosity objects than in our quasar sample, further 
evidence of a radio luminosity and/or redshift dependence 
of clustering strength. Also pertine nt to this is the 2-d egree 
Field (2dF) QSO Redshift Survey ( |Croom et al. 199^ ) cur- 
rently underway at the Anglo- Australian Telescope, which 
will, in due course, measure the clustering of high-z radio- 
quiet quasars directly, and also determine the evolution in 
the clustering strength with redshift. 



8.3 The link between cluster richness and radio 
luminosity 

Does our discovery of a correlation between radio luminosity 
and _B gq imply that environment is the primary factor in 
controlling radio luminosity? If so there are at least three 
ways in which such a situation could come about: (1) the 
environment determines the bulk kinetic power in the radio 
jets, Q, through a correlation of cluster richness with the 
black-hole mass in the central group/cluster galaxy; (2) the 
environment determines Q through more AGN fuel being 
available in richer environments, or (3) Q is independent of 
environment and similar for all radio sources, with the radio 
luminosity being determined solely by the density into which 
the source expands. 

The argument in favour of possibility (1) is based on 
the correlation between black hole mass and the mass of the 
s phero idal component sugge sted b y Kormendy & Richstone 
( 1995| ) and Magorrian et al. (1995). This implies that the gi- 



ant elliptical hosts of radio galaxies and quasars should have 
high black hole masses, ~ 1O 8 -1O 9 M0. If the jets are pow- 
ered by accretion then it seems reasonable to expect that 
the maximum accretion rate is proportional to the Edding- 
ton rate, and hence the black hole mass. So more massive 
galaxies should power more luminous radio sources. These 
massive galaxies will always prefer richer environments, and 
the correlation between luminosity and cluster richness may 
just reflect an increasing mass of the host galaxy (and hence 
of the black hole). The principal argument against this is 
the lack of any obvious correlation of host galaxy luminos- 
ity with ra dio luminosity, at least out to z ~ 0.6 ( Serjeant 
et al. 1997), although there is a suggestion tha t by 2 ~ 1 the 



hosts of radio gal axies do show such a trend ( Roche, Eales 
fc Rawlings 1998| ). 

The possibility (2), that fuelling of the AGN deter- 
mine s its environment is discussed by Ellingson, Green & 
Yee (1991). A group or poor cluster environment may be 
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ideal for the fuelling of a black hole, as encounters will be 
more common than in the field, but will be of low enough 
relative velocity to disrupt the interacti ng ga laxies and cause 
gas to flow into the centre. Lake et al. (1998), however, have 
also demonstated that galaxy harassment (frequent high- 
speed galaxy encounters) in relatively high velocity disper- 
sion clusters may be an efficient mechanism for transporting 
gas to the centre of galaxies. If this mechanism is responsi- 
ble for fuelling the black hole it may explain the correlation 
of radio luminosity with cluster richness, but cannot explain 
why low-power FRI radio sources are also found in rich clus- 
ters. 

Possibility (3) is suggested by models of radio sources 
which show that the radio luminosity of an expanding radio 
source is primarily determined by the bulk kinetic power in 
the radio jet, Q, and the environmental density, p, as a func- 
tion of distance from the nucleus, r, usually parameterised 
by 

p oc po (ao/r) 13 , 

where ao and po are the characteristic radius and den- 
sity, and (5 is a power-law slope. Kaiser, Dennett-Thorpe 
& Alexander (1997) show how the track of the expanding 
radio source on the radio-luminosity - size diagram is de- 
termined by different values of these parameters. The small 
dispersion in galaxy mass (and therefore perhaps black hole 
mass) of radio source host galaxies might suggest that, if all 
radio sources are accreting matter at the same fraction of 
the Eddington rate, the values of Q should be similar for 
all radio sources. However, it is clear from the models that 
radio luminosity is also strongly dependent on the mean en- 
vironmental density. This is because objects located in dense 
gaseous haloes have lower losses due to adiabatic expansion, 
and the transfer of AGN power into radiation thus becomes 
very efficient due t o the reduced expansion losses. Indeed, 
Barthel & Arnaud (1996) use the case of Cygnus A, which 
is in such an environment, to argue that environmental den- 
sity controls radi o sou rce luminosity. Consistent with this, 
Hutchings et al. (1996) claim to find a weak, positive corre- 
lation between lobe luminosity of quasars and B sq . 

On the other hand, a relation between radio and [OIII] 
emission-line luminosity is observed (Rawlings & Saunders 
1991; Willott et al. 1999). The [OIII] A5007 luminosity is be- 
lieved to be powered by UV radiation from the nucleus (for 
a different view, see Dopita & Sutherland 1995), and these 
observations thus argue that radio luminosity is prim arily 
determined by jet kinetic power. Serjeant et al. (1998) also 
find that the radio and optical luminosity of quasars corre- 
late. If both are related to the luminosity of the AGN pow- 
ered by accretion, then this also suggests that jet power, 
not environmental density, is the most important factor in 
determining the radio luminosity. 

To investigate the possible implications of (3) for our 
study though, we have constructed a simple model of how 
radio luminosity should relate to -B gq . We assumed that 
Bgq = c log p and calibrated this relation using the gas den- 
sities from X-ray observations of 21 Abell cluste rs in Jones 
& Forman (|l984|) for which Yee & Lopez-Cruz (|l999|) have 



Including CSS sources 



calculated B cg . The median density and cluster richness for 
these 21 clusters was found to be 3.03xl0 -3 cm -3 and 1075 
Mpc 1 ' 77 . We took this combination of p and B cg to be typ- 
ical for rich clusters. For galaxies with no enhancement in 
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Figure 9. The solid line corresponds to the predicted relationship 
between quasar radio luminosity and environment from simple 
models assuming that environment is entirely responsible for de- 
termining the radio luminosity. The predicted correlation is much 
stronger than what is observed, shown by the dotted line. In this 
plot, we have assumed that logp = B gq /c, where c = 520.6. The 
upper plot shows the relationship when all quasars are included 
in the fit with a radio luminosity-size slope of a = —0.408. In the 
lower plot, the CSS sources in the sample have been ignored, and 
a = —0.239. It may be noted that this rather simple model goes 
very quickly to negative B gq values. This is implicit in the model 
due to the assumption that B gq is proportional to log p. 



t he surrounding galaxy density, we used B gg = 60 Mpc 



(Davis & Peebles 1983) and a density of 3.4 x 10 cm 



which is a typical gas density at 100 kpc in isolated ellip- 
ticals (Forman, Jones & Tucker 1985). This gave a slope of 



, gq -log p relation. 



c = 520.6 for the B g 

We then assumed that L408MHZ oc r a p' 1 ri , as predicted 
by simple analytical models of radio sources (e.g. Miller, 
Rawlings & Saunders 1993). The slope a, was found by fit- 
ting a straight line to the data in a log L408MHz-log r plot, 
which came out as a = —0.408. For consistency, we also 
found the best-fitting slope for the extended (> 30 kpc) 
sources only, which gave a — —0.239. It is interesting to 
note that the quasars show an anticorrelation between radio 
luminosity and size, but that the population of radio sources 
as a whole does not. A likely explanation could be the optical 
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selection effects in the 7C sample. Also the emission at lower 
frequencies might be absorbed in the CSS sources, which 
would explain why there are fewer CSS sources among the 
7C quasars (selected at 151 MHz) as opposed to the MRC 
quasars (selected at 408 MHz) . The correlation between ra- 
dio luminosity and size also weakens when we exclude the 
CSS sources. 

By assuming a fixed jet power, Q, for all the sources, we 
then evaluated the relationship between -£/408MHz and B gq as 

12 

B m / c = — (log L 4 osMHz - a log r) , 

in the two cases of a. The relation is plotted in Fig. ^, where 
the predicted slope of 12/7 is seen to be steeper than what 
the data indicates. So the correlation of _B gq with environ- 
ment is weaker than we expect from simple models in which 
the density alone determines the radio luminosity. We can 
therefore rule out as strong a B gq -Z/408MHz dependence as 
we would see if all objects had about the same jet power, and 
environment was entirely responsible for determining radio 
luminosity. Nevertheless, as the luminosity function for the 
radio jet power is likely to be steeply declining at high pow- 
ers, it seems not unlikely that selection effects could operate 
to produce some correlation between B gq -L408MHz without 
it being as strong as it would be in this rather extreme model 
in which the jet power is the same for all radio sources. 



Given the large scatter in both the B gq -Z/4 



and 



the I/408MHz-A/b (or Lroni]) correlations it is quite possible 
that both environment and radio jet power play important 
roles in determining the radio luminosity. The relationship 
between radio sources and their environments must be com- 
plex, and the vast majority of radio sources may lie in some 
sort of cluster-like environment, from groups of only a few 
galaxies to clusters as rich as Abell class 1 or more. The 
weak correlation we see can probably be explained if we 
assume that radio sources occupy a range of environments 
from fairly isolated ellipticals to rich clusters, and their ra- 
dio luminosities are determined both by the environment 
and the kinetic power in the jets. 



9 SUMMARY 

We have investigated the galaxy environments around 21 
radio-loud, steep-spectrum quasars in the luminosity range 
23.8 < log (L 4 08MHz/WHz" 1 sr _1 ) < 26.7 and redshift range 
0.5 < z < 0.82, and find an excess of galaxies within a 0.5 
Mpc radius centered on the quasar. The galaxy excess in 
each quasar field was quantified using the galaxy-quasar 
cross-correlation amplitude, B gq , and an Abell type mea- 
surement, No. 5. Our estimates of B gq is seen to fall approx- 
imately in the same range as that found for quasar envi- 
ronments by other investigators. We also find that No. 5, al- 
t hough originally intended for fields around radio galaxies 
( [Hill fc Lilly 199l| ), is applicable to quasars fields too and 
that ATq.s for our quasar fields is consistent with No. 5 fo r 



fields around radio galaxies studied by Hill & Lilly (1991). 

The average B gq of our quasar sample was found to be 
265±74 Mpc 1 ' 77 , or in terms of ro, the correlation length 
of the galaxy-quasar cross-correlation function, ro = 23.4 
Mpc. If the galaxies around the quasars were distributed 

60 Mpc 1 ' 77 



and ro « 10 Mpc as found for the galaxy-galaxy correlation 
function in the field (e.g. Davis & Peebles 1983). This shows 
that the quasars are sited in richer than average galaxy en- 
vironments, and two of the quasar fields have B gq above 
740 Mpc 1,77 , considered here to be characteristic for clus- 
ters with Abell richness class 1 or greater. 

Our sample extends previous quasar samples to higher 
redshifts whilst maintaining the same range in radio lu- 
minosity. We compared with _B gq measurements in fields 
around radio-loud quasars from these samples (YG87; 
EYG91; Yee & Ellingson 1993), and detected a weak, but 
significant, correlation between radio luminosity and envi- 
ronment richness. The correlation between B gq and environ- 
ment was found to be much weaker than expected from sim- 
ple models in which gas density alone determines the radio 
luminosity. Probably, the correlations between B gq and envi- 
ronment, and between jet power, Q, and quasar UV/optical 
luminosity can be explained if both Q and environment play 
significant roles in determining the luminosity of a radio 
source. 

Studies of radio galaxies (e.g. Hill & Lilly 1991) have 
found that environmental richness increases with redshift 
and that there is no correlation with radio luminosity, con- 
trary to our observations of radio-loud quasars. The appar- 
ent difference between radio-loud quasars and galaxies would 
cause problems for the unified scheme of radio-loud AGN un- 
less the difference can be explained by selection effects o r a 
mix of two different types of objects (Laing et al. 1994) in 
radio galaxy samples. 

Our main conclusion is that the radio-loud steep- 
spectrum quasars in our sample on average occupy environ- 
ments typical of poor clusters with Abell richness class 0, but 
that quasar environments may cover a wide range in rich- 
ness, from groups of galaxies and poor clusters to clusters as 
rich as Abell class 1 or more. In a conventional dark matter 
cosmology, very rich clusters at redshifts 0.5-0.8 should be 
rare as they would not have had time to form. Nevertheless, 
a few very luminous z ~ 0.8 — 1.3 X-ray selected clusters 
have been found (e.g. Henry et al. 1997; Rosati et al. 1999), 
so at least some very gas-rich environments exist at high 
redshift. It is thus perhaps surprising that so few radio- loud 
objects are found in clusters of Abell richness 1 or higher, 
given the arguments about both galaxy-black hole mass and 
the relationship between gas density and radio luminosity. 
This adds weight to the suggestion that fuelling problems 
may be the reason why radio sources seem to prefer only 
moderately rich environments. 
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